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Review on the inclusive rare decays B Xgj and B Xdj in 
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Abstract. We review the NLL QCD calculations for the branching ratio of _B Xs"/ in the SM. In 
particular, we emphasize the problem related to the definition of the charm quark mass which leads to 
a rather large uncertainty of the NLL predictions. The various steps needed for a NNLL calculation, in 
which the rric issue can be settled, is also sketched. We briefly summarize the results of a calculation of 
the 0{a'^nf) corrections to BR(i3 — > Xsj), which was recently performed as a first step in the NNLL 
program. We then also briefly review the status of the photon energy spectrum and show the comparison 
with experimental data. Finally, we review the status of the CKM suppressed decay mode B Xd'y- 



1 Introduction 



In the Standard model (SM), rare B decays like B — ^ 
or B ^ Xsi^ i~ are induced by one-loop diagrams, where 
virtual W bosons and up-type quarks are exchanged. In 
many extensions of the SM, there are additional contri- 
butions, where the SM particles in the loop are replaced 
by nonstandard ones, like charged Higgs bosons, gluinos, 
charginos etc. If the masses of these new particles are not 
heavier by many orders of magnitude than the heaviest 
SM particles, the new physics contributions to rare B me- 
son decays are expected to be generically large. The sen- 
sitivity for nonstandard effects implies the possibility for 
an indirect observation of new physics, or allows to put 
limits on the masses and coupling parameters of the new 
particles. 

It is obvious that it is only possible to fully exploit 
the new physics potential of these decays when both, pre- 
cise measurements and precise theoretical SM calculations 
exist. 

In the following we mainly concentrate on the decays 
B Xs^ and B X^'j, while the rare semileptonic de- 
cay B Xsi'^i^ is reviewed at this conference by T. 
Hurth [1]. There are experimental analyses of the branch- 
ing ratio BR{B X^j) by CLEO [2-4], ALEPH [5], 
BELLE [6], and BABAR [7] as shown in fig. 1, leading 
to the world average [8] 

BR(B ^ Xs7)oxp = (3.34 ± 0.38) x 10"^ 

In contrast to the exclusive rare decay B K*'-f, 
the inclusive counterpart B — > Xsj is theoretically much 
cleaner as no specific model is needed to describe the 
hadronic final state. Indeed, nonperturbative effects in 
the inclusive decay mode are well under control due to 
the heavy quark expansion technique (HQE), which im- 
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Fig. 1. Branching ratio BR(i3 Xsy): The numbers attached 
to the various experiments reflect the year of publication of 
the corresponding result. The dashed (dotted) band shows the 
theoretical results based on the MS-bar (pole mass) interpre- 
tation of the charm quark mass; see eqs. (4) and (3). Figure 
taken from [7] and world average added from [8]. 



plies that the decay width r{B Xgj) is well approxi- 
mated by the partonic decay rate r{b Xgj) which can 
be analyzed in renormalization group improved pertur- 
bation theory. The (nonperturbative) power corrections 
which scale like l/m^ [9] and l/m^ [10] were estimated to 
be well below 10%. 



2 Theoretical framework 

Short distance QCD effects enhance the partonic decay 
rate r{b — *■ 57) by more than a factor of two. Analytically, 
these QCD corrections contain large logarithms of the 
form a"(mb) In™ (m^/Af), where M = mt or M — mw 
and m < n (with n = 0, 1, 2, ...). In order to get a reason- 
able prediction for the decay rate, it turns out that one 
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has to resum both, the leading-log (LL) terms (m = n) as 
well as the next-to-leading- log (NLL) terms (m = n — 1). 

To achieve the necessary resummations, one usually 
constructs in a first step an effective low-energy theory 
and then resums the large logarithms by renomalization 
group techniques. The low energy theory is obtained by 
integrating out the heavy particles which in the SM are 
the top quark and the VF-boson. The resulting effective 
Hamiltonian relevant for 6 — > 57 in the SM and many of 
its extensions reads 

H,s{b^s^) = -^\tJ2Ci{iJ,)Oi{f,) , (1) 

where Oi(/x) are local operators consisting of hght fields, 
Ct(/u) are the corresponding Wilson coefficients, which 
contain the complete top- and W- mass dependence, and 
At = VtbV*g with Vij being the CKM matrix elements. The 
CKM dependence globally factorizes, because we work in 
the approximation = 0. 

In the basis introduced by Misiak [11], the operators 
read 

^7 = -r:r^'^b{p)s cr'^'' Rb Ff_i^ , 

08 = Y^m(,(M)sa^''i?T»6G«, . (2) 

As the Wilson coefficients of the QCD penguin operators 
O3, ...,06 arc small, we do not list them here. 

A consistent calculation for 6 — > 57 at NLL precision 
requires three steps: 

1) a matching calculation of the full standard model the- 
ory with the effective theory at the scale /x = to 
order al for the Wilson coefficients, where iJ,w denotes 
a scale of order Myy or rrit; 

2) a renormalization group evolution of the Wilson coef- 
ficients from the matching scale jiw down to the low 
scale /Ufe = 0{rrH)), using the anomalous-dimension ma- 
trix to order a^; 

3) & calculation of the matrix elements of the operators 
at the scale /x = /Xb to order a\ . 

As all three steps are rather involved, a common ef- 
fort of several independent groups was needed in order to 
calculate the NLL prediction for BR(S Xsj) [12 15, 
11,16-18]. For a detailed summary of the various steps 
and intermediate results, we refer to the recent review by 
T. Hurth [19]. However, we would like to point out that 
the most difficult part, viz. the calculation of three- loop 
anomalous dimensions performed by Chetyrkin, Misiak 
and Miinz in 1996 [11], was only confirmed very recently 
by Gambino, Gorbahn and Haisch [20]. Their paper also 
contains the three- loop mixing of the four- Fermi operators 
into O9, which is important for the process B Xs£~^i~. 

During the completion of the NLL QCD corrections, 
also calculations of electroweak corrections were started 
[21-23]. At present, the corrections of order agm ln(/x&/M) 
[ag ln(/i5/M)]", as well as the sublcading terms of order 
ctein[cts^^{l^b/M)]"- [24] are Systematically available. 



3 NLL (and partial NNLL) results for 
BR{B X«7) 

Combining NLL QCD corrections with the electroweak 
corrections just mentioned and also including the 1/m^ 
[9] and l/m^ [10] power corrections, the branching ratio 
reads 

BR{B Xs7) = (3.32 ± 0.14 ± 0.26) x lO"'' , (3) 

where the first error reflects the dependence on the renor- 
malization scale /ib varied in the interval mb/2 < /it, < 
2mb, while the second error reflects the error due to the 
uncertainties in the input parameters. 

Among the input parameters the charm quark mass 
rric plays a crucial role. The charm quark mass depen- 
dence only enters the prediction for the decay width at the 
NLL level, more precisely through the 0{as) correction to 
the matrix elements {sj\Oi^2\b) ■ Until recently, all authors 
used the pole mass value mj?"''' for the charm quark mass 
in numerical evaluations, leading to a branching ratio as 
specified in eq. (3). 

In 2001, however, Gambino and Misiak [25] pointed out 
that the MS-bar mass ttTc, normalized at /i « mb/2, could 
be the better choice, because the charm quark appears as 
an off-shell particle in the loop involved in the above men- 
tioned matrix element with a typical virtuality of mb/2. 
Using this interpretation, rn.c/mi, = 0.22 ±0.04 is substan- 
tially smaller than the value m^°^'^/mi, = 0.29 ± 0.02 used 
in eq. (3), leading to a branching ratio of [25,26] 

BR(B ^ X,7) = (3.70 ± 0.30) x 10"^ . (4) 

We would like to stress here that the above argument 
in favour of Trie is an intuitive one. Formally, the difference 
between using mP°'® or mc amounts to a NNLL effect at 
the level of the branching ratio. This means that a NNLL 
becomes necessary in order to unambiguously fix this is- 
sue. 

Before sketching the NNLL program, we would like to 
stress that settling the issue is also important when 
extracting bounds on new physics, based on NLL cal- 
culations [27,28]. For example, in the type-II two-Higgs- 
doublet model, one obtains a bound from 6 — > S7 on the 
charged Higgs boson mass of niH > 350 GeV (99% C.L.) 
when using rrTc- When using on the other hand m^°^^, the 
bound is TTiH > 280 GeV (99% C.L.) [25]. 

Concerning the NNLL program, it is clear that in or- 
der to get a full NNLL QCD result for BR{B -> Xsj), 
all the three steps listed above have to be improved by 
one order in . This means that three-loop matching cal- 
culations are needed, up to four-loop anomalous dimen- 
sions have to be worked out and up to three-loop calcu- 
lations at the level of the matrix elements {s"f\Oi{iJ,b)\b) 
have to be performed. Several groups have been formed 
in order to attack this ambitious goal. Recently, a cal- 
culation of the 0{agnf) corrections to the matrix ele- 
ments of the operators Oi , O2, O7 and Og was published 
[29]. Diagrammatically, these contributions are generated 
by inserting quark bubbles (n/ denotes the number of 



K. Bieri and C. Greub: B X^^ and B Xd7 in the SM 



3 




m 1.5 

1 

0.5 

\ 

3 4 5 6 7 8 9 

H (GeV) 

Fig. 2. BR(i3 Xs"t) as a function of the renorniahzation 
scale 11. The dash-dotted curve shows the NLL prediction; the 
solid curve contains in addition the ©(a^n/) terms. The long- 
dashed (short-dashed) curve is obtained by switching off the 
C(ajn/) corrections to O7 (O2). Figure taken from ref. [29]. 



light quarks) into the gluon propagators in the diagrams 
which are involved in the calculations at NLL order. We 
note that these contributions are not related to the defi- 
nition problem of mg. However, in many other cases they 
are sources of large corrections. E.g., in the scmilcptonic 
decay width r{B Xcii^g) these ©(a^n/) terms (af- 
ter replacing nj ^ — 3/3o/2, according to the procedure 
of naive non-abelianization) incorporate more than 80% 
of the complete 0(a.'^^) corrections [30]. The impact of 
the 0{a^nf) corrections to BR(i? ^.sT) ^u'^' shown in 
fig. 2. The dash-dotted curve shows the NLL prediction, 
while the solid curve incorporates in addition the 0{a'lnf ) 
terms (after the replacement Uf ^ — 3/3o/2, according to 
naive non-abelianization). As one sees from the figure, the 
0{a'^nf) corrections seem to be small. Note, however, that 
this is a result of a relatively large accidental cancellation 
between corrections to O2 and O7. This point is illustrated 
by the long-dashed and short-dashed curve, which are ob- 
tained by switching off the (^(a^n/) corrections to O7 and 
O2, respectively. 



4 Partially integrated BR and photon energy 
spectrum 

The photon energy spectrum of the partonic decay b sj 
is a delta function, concentrated at ~ {nib/ 2), when the 
6-quark decays at rest. This delta function gets smeared 
when considering the inclusive photon energy spectrum 
from a B meson decay. There is a pcrturbative contri- 
bution to this smearing, induced by the Bremsstrahlung 
process b — > s^g [16, 17], as well as a nonperturbative one, 
which is due to the Fermi motion of the decaying b quark 
in the B meson. 

For small photon energies, the 7-spectrum from B — > 
Xs"f is completely overshadowed by background processes, 
like b cudj and b uud'j. This background falls off 
very rapidly with increasing photon energy, and becomes 
small for > 2 GeV [31]. This implies that only the 
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Fig. 3. Photon energy spectrum: The data points represent 
the CLEO measurement [4]. The histogram shows the theory 
result based on the spector model using pp = 410 MeV and 
(mb) = 4.69 GeV. Fig. taken from ref. [4]. 

partial branching ratio 

BR(B-X,7)«")=y^_^^_^ -^dE, (5) 

can be directly measured, with E^^^ = 0(2) GeV. 

Putting the energy cut at = 2.0 GeV, CLEO 

used two methods to analyze their data on the photon en- 
ergy spectrum in their most recent analysis: First, the Ali- 
Grcub model [16,32], based on the spectator model for- 
mulated in ref. [33] and second, methods based on HQET 
[22]. 

The spectator model contains two free parameters, viz. 
Pf, the average Fermi momentum of the b quark in the B 
meson and the mass of the spectator quark, rrispec- Equiva- 
lently {pp, (rrib)) can be used as the free parameters, where 
(rrib) is the average b quark mass as defined in ref. [16,32]. 
In Fig. 3 a comparison between theory and experiment is 
shown. Using pp = 410 MeV and (mb) = 4.69 GeV the 
best fit is obtained. We would like to stress that similar 
values for these parameters are also obtained when fitting 
the lepton spectra in _B ^ Xc£i^ and B Xu^v- 

A modern way - based on first principles - implements 
the Fermi motion in the framework of the heavy-quark 
expansion. When probing the spectrum closer to the end- 
point, the OPE breaks down, and the leading twist non- 
perturbative corrections must be resummed into the B 
meson structure function /(fc+) [34], where fc+ is the light- 
cone momentum of the b quark in the B meson. The phys- 
ical spectrum is then obtained by the convolution 

dk^f{K)^{ml) , (6) 

where {drpg,rt/dEy){ml) is the partonic differential rate, 

written as a function of the "effective mass" ml = ?Tif,-|-fc+. 
The function f{k^) has support in the range —00 < A:+ < 
A, where A = tub — rub in the infinite mass limit. This im- 
plies that the addition of the structure function moves the 
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Fig. 4. Partially integrated branching ratio as a function of 
the energy cutoff iS™'"; Curves taken from Kagan and Neubert 
[22] . Data point represent CLEO measurements [2-4] . 



not small relative to At and Ac; therefore, also the current- 
current operators O" and O2, weighted by A„, contribute. 




Fig. 5. Interaction of the c— and m— quark loop with soft glu- 
ons surrounding the b quark in the B meson. 



partonic cndpoint of the spectrum from TOfc/2 to the physi- 
cal endpoint 1113/2. While the shape of the function f{k+) 
is unknown, the first few moments An ~ J dk^ k^f{k+) 
are known: Aq = 1, Ai = and A2 = — Ai/3. As An 
(n > 2) are poorly known, several Ansatze were used 
for f{k+); e.g. Neubert and Kagan [22] used f{k+) = 
N{1 - a;)''e(^+'')^, with x = k+/A. Taking into account 
the constraints from ^1 and A2, the independent pa- 
rameters in this Ansatz can be chosen to be 777,5 and Ai. 
As shown in [22], the uncertainty of m,b dominates the er- 
ror of the partial branching ratio. In figure 4 the partial 
branching ratio is shown for the_ relevant range of mb as 
a function of E™^'^, keeping Xi/A"^ fixed. The data points 
show three CLEO measurements. In the oldest one the 
photon energy cut was put at £^™'" = 2.2 GeV, while in 
the most recent analysis this cut was lowered to 2.0 GeV, 
which is very important, because at 2.0 GeV the theoret- 
ical error on the partial branching ratio is considerably 
smaller, as seen from fig. 4. 

To determine from the measurement of the partial 
branching ratio the full BR, one needs from theory the 
fraction R of the B Xs"f events with photon energies 
above E^'". Based on Kagan-Neubert [22], CLEO [4] ob- 
tained R = (O.QlSlgQgg) . A similar result is also obtained 
when using the spectator model. 

It has been shown that up to corrections of 0{Aqcd / 'nib] 
the same shape function also describes B — > X^iv [35]. 
This implies that the photon energy spectrum can be 
used to predict the fraction oi B ^ Xuiv events with 
iJiept > 2.2 GeV, where leptons coming from B — > X^tv 
are absent for kinematical reasons. Taking into account 
perturbative and ^Iqcd /to6 corrections [36-38] , it is pos- 
sible to extract Vuh from a measurement of the B Xuf-iy 
decay rate in the region above 2.2 GeV. CLEO used this 
strategy in ref. [39] to extract the CKM matrix element 
\Vub\, obtaining \Vub\ = (4.08 ± 0.56exp ± 0.29th) x lO'^. 



5 .B ^ Xd-y in the SM 

The decay B Xd"f can be treated in a similar way as 
B Xs'y [40]. The only difference is that A„ for 6 ^ ^7 is 



Unfortunately, these operators induce long-distance con- 
tributions to B ^ Xd'y, which at present are not very well 
understood. To illustrate the problem, we first look at the 
corresponding charm quark loop, depicted in fig. 5. In this 
case, one can expand the loop function 



dx 



dy- 



xy 



l-^x{l-x)-2xy^-^ 



in powers of t = kgkj/m'^, where kg and k-y denote the 
momentum of the gluon and the photon, respectively. This 
expansion generates the so-called Voloshin terms [10] . which 
in BR(i3 Xgj) are a 3% effect. Obviously, there is 
no such OPE in the case of the u— quark loop. However, 
Buchalla, Isidori and Key [41] argued that an expansion 
in l/t can be done, leading to non-local operators. From 
naive dimensional counting, the leading contribution is ex- 
pected to be of order ylgcD/^^h- 

In reference [40], where NLL calculations for the pro- 
cess B Xd"f were presented, the uncertainties due to 
the long-distance effects were absorbed into the theoret- 
ical error. Using Hb = 2.5 GeV and the central values 
of the input parameters, the analysis in reference [40] 
yields a difference between the LL and NLL predictions 
for BR(i? Xd'j) of ~ 10%, increasing the branching 
ratio in the NLL case. For a fixed value of the CKM- 
' Wolfenstein parameters p and rj, the theoretical uncer- 
tainty of the average branching ratio {BR{B Xdj)) of 
the decay B Xd'y and its charge conjugate B — s- Xd'y 
is: A{BR{B Xd7))/(BR(B ^ X^i)) = ±(6- 10)%. Of 
particular theoretical interest for constraining p and r/ is 
the ratio of the branching ratios, defined as 



Ridj/sj) 



{BR{B ^ Xdj)) 



(7) 



(BR(S ^ Xs-f)) ' 

in which a good part of the theoretical uncertainties can- 
cels. Varying the CKM- Wolfenstein parameters p and rj in 
the range -0.1 < p < 0.4 and 0.2 < ?? < 0.46 and tak- 
ing into account other parametric dependences, the results 
(without electroweak corrections) are 

6.0 X 10"*^ < BR{B Xd-f) < 2.6 x 10"^ , 
0.017 < Ridj/sj) < 0.074 . 
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Another observable, which is also sensitive to the CKM 
parameters p and rj, is the CP rate asymmetry acp, de- 
fined as 

^ ^ r{B^Xaj)-r(B^xh) 

r{B ^ Xrf7) + r{B ^ Xd7) ■ ^ 

Varying p and ry in the range specified above, one gets 
7% < acp < 35% [40]. We would like to point out that 
acp is at the moment only available to LL precision and 
therefore suffers from a relatively large renormalization 
scale dependence. 

In summary, this decay mode is very challenging, both 
in theory and experiment: On the theory side more work is 
needed concerning the nonperturbative contributions as- 
sociated with the u— quark loop, while on the experimental 
side the observation of this decay needs high statistics and 
a very good discrimination between pions and kaons. 

This work is partially supported by the Swiss National Founda- 
tion and by RTN, BBW-Contract NO. 01.0357 and EC-Contract 
HPRN-CT-2002-00311 (EURIDICE). 
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